Partial dynamical symmetry in quantum Hamiltonians with higher-order terms 
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A generic procedure is proposed to construct many-body quantum Hamiltonians with partial 
dynamical symmetry. It is based on a tensor decomposition of the Hamiltonian and allows the 
construction of a hierarchy of interactions that have selected classes of solvable states. The method 
is illustrated in the SO (6) limit of the interacting boson model of atomic nuclei and applied to the 
nucleus 196 Pt. 
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The concept of dynamical symmetry (DS) is now 
widely accepted to be of central importance in our un- 
derstanding of many-body systems. In particular, it 
had a major impact on developments in nuclear [l| and 
molecular 0] physics and made significant contributions 
to virtually all areas of many-body physics Its ba- 
sic paradigm is to write the Hamiltonian of the system 
under consideration in terms of Casimir operators of a 
set of nested algebras. Its hallmarks are (i) solvability 
of the complete spectrum, (ii) existence of exact quan- 
tum numbers for all eigenstates, and (iii) pre-determined 
symmetry-based structure of the eigenfunctions, inde- 
pendent of the Hamiltonian's parameters. 

The merits of a DS are self-evident. However, in most 
applications to realistic systems, the predictions of an ex- 
act DS are rarely fulfilled and one is compelled to break 
it. More often one finds that the assumed symmetry is 
not obeyed uniformly, i.e., is fulfilled by only some states 
but not by others. The need to address such situations 
has led to the introduction of partial dynamical symme- 
tries (PDSs). The essential idea is to relax the stringent 
conditions of complete solvability so that the properties 
(i)-(iii) are only partially satisfied. Partiality comes in 
three different guises: (a) part of the eigenspectrum re- 
tains all the DS quantum numbers 0, (b) the en- 
tire eigenspectrum retains part of the DS quantum num- 
bers [y, 0] , and (c) part of the eigenspectrum retains part 
of the DS quantum numbers [8|. PDS of various types 
have been shown to be relevant to nuclear [j| |(| 0, M, @ 
and molecular [Io| spectroscopy, to systems with mixed 
chaotic and reg ular dynamics and to quantum phase 
transitions [121 ]. 

We emphasize that DS as well as its generalization PDS 
are notions that are not restricted to a specific model but 
can be applied to any quantal system of interacting par- 
ticles, bosons and fermions. In this Letter we propose a 
generic method to construct quantum Hamiltonians with 
PDS of type (a) and show that its existence is closely re- 
lated to the order of the interaction among the particles. 
The procedure is discussed in general terms and subse- 
quently illustrated with a nuclear-physics example with 
an application to the nucleus 196 Pt. 



The analysis starts from the chain of nested algebras 
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where, below each algebra, its associated labels of irre- 
ducible representations (irreps) are given. Eq. ([1]) im- 
plies that Gdyn is the dynamical algebra of the system 
such that operators of all physical observables can be 
written in terms of its generators [l3l |; a single irrep of 
Gdyn contains all states of relevance in the problem. In 
contrast, G sym is the symmetry algebra and a single of 
its irreps contains states that are degenerate in energy. 
A frequently encountered example is G sym = SO (3), the 
algebra of rotations in 3 dimensions, with its associated 
quantum number of total angular momentum L. Other 
examples of conserved quantum numbers can be the total 
spin S in atoms or total isospin T in atomic nuclei. 

The classification ((T|) is generally valid and does not 
require conservation of particle number. Although the 
extension from DS to PDS can be formulated under such 
general conditions, let us for simplicity of notation as- 
sume in the following that particle number is conserved. 
All states, and hence the representation [h], can then be 
assigned a definite particle number N. For N identical 
particles the representation [h] of the dynamical algebra 
Gdyn is either symmetric [N] (bosons) or antisymmetric 
[1^] (fermions) and will be denoted, in both cases, as 
[ft. at]. For particles that are non-identical under a given 
dynamical algebra Gd yn , a larger algebra can be chosen 
such that they become identical under this larger algebra 
(generalized Pauli principle) [l4j . The occurrence of a DS 
of the type dJ) signifies that eigenstates can be labeled 
as | [h n] (E) ... A); additional labels (indicated by . . .) are 
suppressed in the following. Likewise, operators can be 
classified according to their tensor character under ([1]) 
as ]( CT ) A . Of specific interest in the construction of a 
PDS associated with the reduction |T]), are the n-particle 
annihilation operators T which satisfy the property 



f {hn]{a)x \[h N ](Z )A) = 0, 



(2) 



for all possible values of A contained in a given irrep (Eo). 
Any number-conserving normal-ordered interaction writ- 



2 



ten in terms of these annihilation operators (and their 
Hermitian conjugates which transform as the correspond- 
ing conjugate irreps) can be added to the Hamiltonian 
with a DS (Q]), while still preserving the solvability of 
states with (E) = (So)- The annihilation condition ((2|) 
is satisfied if none of the G irreps (E) contained in the 
Gdyn irrep [/ijv— «] belongs to the G Kronecker product 
(So) x (a). So the problem of finding interactions that 
preserve solvability for part of the states (P) is reduced 
to carrying out a Kronecker product. 

Let us now illustrate this procedure with an exam- 
ple taken from nuclear physics where the concept of DS 
has been applied successfully in the context of the in- 
teracting boson model (IBM) In this model, low- 
energy collective states are described in terms of N in- 
teracting monopole (s) and quadrupole (d) bosons repre- 
senting valence nucleon pairs. The dynamical algebra is 
Gdyn = U(6) and the symmetry algebra is G sym = SO (3). 
Three DS limits occur in the model with leading subal- 
gebras U(5), SU(3), and SO(6), corresponding to typical 
collective spectra observed in nuclei, vibrational, rota- 
tional, and 7-unstable, respectively. Here we focus on 
the SO (6) limit with predictions that were found to cor- 
respond closely to the empirical structure of some plat- 
inum nuclei [T5| as well as of other nuclei, notably, around 
mass number A = 130 [l6j . 

The classification of states in the SO(6) limit of the 
IBM is 

U(6) D SO(6) D SO(5) D SO(3) D SO(2) 

I I I i I ■ (3) 

[N] (E) (r) v A L M 

The multiplicity label v A in the SO (5) D SO (3) reduction 
will be omitted in the following when it is not needed. 
The eigenstates \[N](Yi)(t)i/aLM) are obtained with a 
Hamiltonian with SO (6) DS which, for one- and two- 
body interactions, can be transcribed in the form 

H US = AP + P-+BC so{5) +CC so{3) . (4) 

Here Cq denotes the quadratic Casimir operator of G, 
P + = I( s t s t _ rft . rft), 4P+A. = N(N + 4) - G SO ( 6 ) 

and P- — PI. The total boson number operator, N = 
n s +nd, is the linear Casimir of U(6) and is a constant for 
all iV-boson states. The spectrum of i?DS is completely 
solvable with eigenenergies 

E DS = iA(iV-E)(iV + E + 4) + Sr(T + 3) 

+ CL(L + 1). (5) 

In contrast, Hamiltonians with SO(6) PDS preserve 
the analyticity of only a subset of the states ([3]). The 
construction of interactions with this property requires 
n-boson creation and annihilation operators with definite 
tensor character in the basis ([3]): 

B[ n ](<7)(T)lm> ^[n 5 ](<T>(r)im = ™ (^[n]( ( j)(r);,-m) ' 

(6) 



TABLE I: Normalized two- and three-boson SO(6) tensors. 
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Of particular interest are tensor operators with cr < u. 
They have the property 

B [n5]{a)(T)lm \[N]{N)(T)y A LM)=0, a < n, (7) 

for all possible values of t, L contained in the SO (6) ir- 
rep (N). This is so because the action of B\ n 5u a \i T \i m 
leads to an (N — n)-boson state that contains the SO (6) 
irreps (E) = (N — n — 2i), i = 0, 1, . . . which cannot 
be coupled with (cr) to yield (E) = (N), since a < n. 
Number-conserving normal-ordered interactions that are 
constructed out of such tensors with a < n (and their 
Hermitian conjugates) thus have | [N] (N){t)v a LM) as 
eigenstates with zero eigenvalue. 

A systematic enumeration of all interactions with this 
property is a simple matter of SO (6) coupling. For one- 
body operators one has 

and no annihilation operator has the property ([7]). 
Coupled two-body operators are of the form 

# cx V VG <,T>(r) ' f6 f x6M (i) (9) 

TkT k i kk' 

where (• x •)„ denotes coupling to angular momentum 
I and the G-coefficients are known SO(6) D SO(5) D 
SO(3) isoscalar factors [Hj]. This leads to the normalized 
two-boson SO(6) tensors shown in Table [IJ There is one 
operator with a < n = 2 and it gives rise to the following 
SO(6)-invariant interaction 

^[2]<0)(0)00^[2 5 ]<0>(0)00 = \P+P-, (10) 
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which is simply the SO (6) term in Pds, Eq. ([4]). 
This proves that a two-body interaction which is di- 
agonal in \[N](N)(t)i/ a LM) is diagonal in all states 
\[N](E)(t)v a LM). This result is valid in the SO(6) limit 
of the IBM, but not in general. For example, from a ten- 
sor decomposition of two-boson operators in SU(3) one 
concludes that the SU(3) limit of the IBM does allow a 
PDS with two-body interactions [HQ. 

Three-body operators with good SO (6) labels can be 
obtained from an expansion similar to © and this leads 
to the normalized three-boson SO (6) tensors shown in 
Table HI In terms of the P + operator introduced above, 
the two operators with a < n = 3 are 

fit — Ip.wt Rt - Ip, ,t fill 

D [3](l)(l)2m — 2 Jr + U m> - D [3](l)(0)00 " 2 r + b ' ^ 1L ) 

and from these one can construct the interactions with 
an SO(6) PDS. The only three-body interactions that 
are partially solvable in SO (6) are thus P+h s P- and 
P+hdP—, involving the s- and c£-boson number operators. 
Since the combination P+(n s + n<j)P_ = (N — 2)P + P_ is 
completely solvable in SO (6), there is only one genuine 
partially solvable three-body interaction which can be 
chosen as P + n s P_, with tensorial components a = 0, 2. 

The generalization to higher orders now suggests itself. 
For example, four-body interactions with SO (6) PDS are 
written in terms of B\ m){r)lm and Pf 4]<0)(0)00 , and Her- 
mitian conjugate operators. Without loss of generality, 
these operators can be written as 

^[4](2)(-r)*m 00 ^+^[2](2)(r)im' ^[4]{0)(0)00 K P +' ( 12 ) 

A four-body interaction with SO (6) PDS is thus of the 
form P+V^P- where Vz is an arbitrary two-body in- 
teraction. This interaction leaves solvable the class of 
states with £ = N but, in general, admixes those with 
£ < N. The conclusion is that we can construct a hier- 
archy of interactions P+n s P^, P+V-jP*, P\_V p P k _ of 
order 2k + 1, 2k + 2, . . . , 2k +p, respectively, that retain 
the SO(6)-DS in selected states. If V p is SO(5)-invariant, 
r is a good quantum number for all states. 

The SO(6)-DS spectrum of Eq. resembles that of 
a 7-unstable deformed rotor, where states are arranged 
in bands with SO(6) quantum number £ = N — 2v, 
(v = 0, 1, 2, . . .). The in-band rotational splitting is gov- 
erned by the SO(5) and SO(3) terms in P DS Q. A com- 
parison with the experimental spectrum and E2 rates of 
196 Pt is shown in Fig. Q] and Table M The SO(6)-DS 
limit is seen to provide a good description for properties 
of states in the ground band (E = N) . This observation 
was the basis of the claim [l5j] that the SO(6)-DS is man- 
ifested empirically in 196 Pt. However, the resulting fit to 
energies of excited bands is quite poor. The 0^~, 0g , and 
0j levels of 196 Pt at excitation energies 0, 1403, 1823 keV, 
respectively, are identified as the bandhead states of the 
ground (v = 0), first- (v = 1) and second- (v = 2) excited 
vibrational bands [la ]. Their empirical anharmonicity, 
defined by the ratio R = E(v = 2) /E(v = 1) -2, is found 



TABLE II: Observed [3 and calculated B(E2) values (in 
e 2 b 2 ) for 196 Pt. For both the exact (DS) and partial (PDS) 
SO (6) dynamical symmetry calculations, the E2 operator is 
e b [(s f x d + d 1 " x s) (2) x d) {2) ] with e b = 0.151 eb and 
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to be R = -0.70. In the SO(6)-DS limit these bandhead 
states have r = L — and £ = iV, iV — 2, — 4, respec- 
tively. The anharmonicity R — — 2/(JV+l), as calculated 
from Eq. ([5|), is fixed by N. For = 6, which is the ap- 
propriate boson number for 196 Pt, the SO(6)-DS value is 
R = —0.29, which is in marked disagreement with the 
empirical value. A detailed study of double-phonon exci- 
tations within the IBM, has concluded that large anhar- 
monicities can be incorporated only by the inclusion of 
at least cubic terms in the Hamiltonian [13]. In the IBM 
there are 17 possible three-body interactions. One is thus 
confronted with the need to select suitable higher-order 
terms that can break the DS in excited bands but pre- 
serve it in the ground band. This can be accomplished 
by the PDS construction presented in this work. On the 
basis of the preceding discussion we propose to use the 
following Hamiltonian with SO(6)-PDS 

PpDS = PDS+^P+n s P_, (13) 

where the terms are defined in Eqs. g]) and dUJ) . The 
spectrum of Ppds is shown in Fig. [TJ The states be- 
longing to the £ = N = 6 multiplet remain solvable 
with energies given by the same DS expression, Eq. (O. 
States with £ < 6 are generally admixed but agree better 
with the data than in the DS calculation. For example, 
the bandhead states of the first- (second-) excited bands 
have the SO(6) decomposition £ = 4: 76.5% (19.6%), 
£ = 2: 16.1% (18.4%), and £ = 0: 7.4% (62.0%). Thus, 
although the ground band is pure, the excited bands ex- 
hibit strong SO(6) breaking. The calculated SO(6)-PDS 
anharmonicity for these bands is R = —0.63, much closer 
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FIG. 1: Observed spectrum of 196 Pt pl| compared with the calculated spectra of Hos Q, with SO(6) dynamical symmetry 
(DS), and of -ffpDS l|13p with partial dynamical symmetry (PDS). The parameters in Hus (Hpus) are A — 174.2 (122.9), 
B = 44.0 (44.0), C = 17.9 (17.9), and r\ = (34.9) keV. The boson number is N = 6 and E is an SO(6) label. 



to the empirical value, R = —0.70. We emphasize that 
not only the energies but also the wave functions of the 
£ = N states remain unchanged when the Hamiltonian is 
generalized from DS to PDS. Consequently, the E2 rates 
for transitions among this class of states are the same in 
the DS and PDS calculations. This is evident in Table HTl 
where most of the E2 data concern transitions between 
X = N == 6 states. Only transitions involving states 
from excited bands (e.g., the Ojj" state in Table [TTJ) can 
distinguish between DS and PDS. Unfortunately, such in- 
terband E2 rates are presently poorly known experimen- 
tally. Their measurement is highly desirable for further 
testing the SO(6)-PDS wave functions. 

In conclusion, we have presented a systematic proce- 
dure for identifying and selecting interactions, of a given 
order, with partial dynamical symmetry (PDS). This al- 



lows the construction of Hamiltonians that break the dy- 
namical symmetry (DS), but retain selected subsets of 
solvable eigenstates with good symmetry. As demon- 
strated in this work, the advantage of using higher-order 
interactions with PDS is that they can be introduced 
without destroying results previously obtained with a DS 
for a segment of the spectrum. These virtues generate an 
efficient tool which can greatly enhance the scope of al- 
gebraic modeling of quantum many-body systems. 
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